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Abstract 

We define, within the Schrodinger functional scheme (SF), the matrix el¬ 
ements of the twist-2 operators corresponding to the first two moments of 
non-singlet parton densities. We perform a lattice one-loop calculation that 
fixes the relation between the SF scheme and other common schemes and 
shows the main source of lattice artefacts. This calculation sets the basis for 
a numerical evaluation of the non-perturbative running of parton densities. 
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1 Introduction 


The evolution with the energy scale of hadronic structure functions can be computed 
within renormalization group improved perturbative QCD, but the calculation of 
their absolute normalization needs non-perturbative techniques. The lattice ap¬ 
proach is a theoretically well established framework for non-perturbative estimates 
and various calculations of the hrst moments of hadronic structure functions have 
appeared in the literature |I[. The results show a discrepancy with experimental 
data already for the non-singlet structure functions |^, There are various sources 
of errors that may explain the disagreement. On the lattice side, possible systematic 
errors are the quenched approximation and the still rather large values of the bare 
coupling constant used in present simulations, or equivalently the rather low values 
of the lattice momentum cutoff, which prevent a correct estimate of the continuum 
limit. On the experimental side, structure functions are reliably parametrized by 
leading-twist operators only at large momentum transfers; a comparison at energy 
scales of a few GeV is affected by systematic errors due to higher-twist effects. A 
standard procedure is then to evolve the lattice results from the cutoff scale, a couple 
of GeV, to higher scales, where the comparison is safe. The evolution is done by us¬ 
ing perturbative results for the renormalization constants of leading-twist operators 
on the lattice, which may be affected by the truncations of the perturbative series. 

Non-perturbative estimates of renormalization constants have already been 
performed for local lowest-dimensional currents [Q and have been seen to be essen¬ 
tial for a proper definition of the quark mass through chiral Ward identities or for 
disentangling the operator mixing entering the evaluation of hadronic matrix ele¬ 
ments of the effective four-fermion weak Hamiltonians |p. A possible method for 
evaluating these renormalization constants is based on the Schrodinger functional 
(SF): this method has been successfully applied to the study of non-perturbative 
running of the coupling constant and of the quark mass 0. We intend to per¬ 
form a similar non perturbative estimate of the running of leading-twist operators 
of hadronic structure functions in the SF scheme. Non-perturbative estimates of 
the running of parton densities are needed if the scale at which the hadron matrix 
element is renormalized is held fixed while the continuum limit is taken. In this 
case the ratio of the renormalization scale over the lattice cutoff becomes large, and 
a perturbative expansion of the leading logarithms of such a ratio is a priori not 
reliable. 

The present paper defines the SF matrix elements of the first and second 
moments of non-singlet leading-twist operators and reports the results on the one- 
loop perturbative calculation. The calculation is useful to learn about the size of 
lattice artefacts and to £x those finite constants peculiar to the SF scheme that can 
be used to make contact with more common schemes such as minimal subtraction. 
In the first section we recall the basic features of the Schrodinger functional involving 
fermions, in the second we discuss the choice of the matrix element for twist-2 non¬ 
singlet operators, and in the third we report the results of our calculations and some 
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concluding remarks. 


2 The Schrodinger functional for fermions 

This section is only meant to recall some basic facts that have been discussed ex¬ 
haustively in the literature [0 . The Schrodinger functional represents the amplitude 
for the time evolution that takes into account quantum fluctuations of a classical 
field configuration between two predetermined classical states. It takes the form of 
a standard functional integral with fixed boundary conditions. It has been shown 
that its renormalizability properties are the same as those of the theory in an infinite 
volume, modulo the possible presence of a finite number of boundary counterterms. 
In QCD with fermions, it can be written as: 


Z\C'J, p'; C, p,p] = j ( 1 ) 

where C\C and p',p',p,p are the boundary values of the gauge and fermion fields 
respectively. In the following discussion the classical boundary gauge held will be 
set to zero. According to ref. @] expectations values may involve the response to a 
variation of the classical Fermi held conhgurations on the boundaries (: 


C(x) - , 

^p(x) 

C(x) = - 

(5 

6p{x.) 

(2) 

6 


6 

(3) 

C (x) - , 

(5p'(x) 

C'(x) = 

V'(x) 


as well as the huctuating Fermi helds between the boundaries tjj. In the following we 
will use the same notations as in ref. ^j. By working out the functional integration 
in the presence of external sources one obtains the expression for the generating 
functional for all possible expectation values involving fermions, and in particular 
the basic non-zero contractions: 


= s{x,y) (4) 

[^(a^)C(y)]F = S{x,y)Uo{y - a6)-^P+\yo=a (5) 

[^(a^)C'(y)]F = S{x,y)Uoiy)P-\yo=T-a (6) 
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[C(x)V’(|/)]f = P-Uo{x - a^)S{x,y%^=a (7) 

[C'(x)^(|/)]f = P+Uo{x)~^S{x,y%^=T-a ( 8 ) 

[C(x)C(y)]F =P-Uo{x - aO)S{x,y)Uo{y - aO)P+\a:o=yo=a 

1 _ (9) 

- 2^-7fc(vI + Vfc)a“^5xy 

[C(x)C'(y)]F = P-Uo{x - aO)S{x,y)Uo{y)P_\^^=a,yo=T-a (10) 

[C'(x)C(y)]F = P+Uo{x)-^S{x,y)Uo{y - a6)-^P+\^,=T-a,yo=a (11) 

[C'(x)C'(y)]F =P+Uo{x)-^S{x,y)Uo{y)P-\^^=y^=T-a 

1 _ (12) 
~ 2'^+7fc(vI “I" Vfc)® '^xy 

where P± = 1(1 ± 70 ) and S{x,y) is the fermion propagator satisfying: 

{D + mo)S{x,y) = a~^5^y, 0 < Xq < T. (13) 


The operator D is the lattice version of the covariant derivative. In this paper 
we will use the Wilson formulation where D takes the form: 

7^ = ^{7/.(vl +Vm)-« vi V/.} (14) 

To one-loop approximation, the above expression must be expanded in pertur¬ 
bation theory up to order g'^ by expanding the gluon unitary matrix U in terms of 
vector potentials q: 

Uy{x) = {1 +gaqy{x)+ 0{g‘^)} ( 15 ) 
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3 The correlations for non-singlet twist-2 opera¬ 
tors 

In the continuum, moments of non singlet structure functions are related, through 
the operator product expansion, to hadronic matrix elements of local twist-2 oper¬ 
ators of the form: 


= {^T V(2:)7[/.i Y'^ix) + trace terms (16) 

<—^ 

where is the covariant derivative, [“indices”] means symmetrization, and p is the 
hadron momentum that can be assumed light-like up to power corrections of order 
with Q the four-momentum transfer of the vector boson probing the hadron 
structure. 

The twist is defined from the difference between the engineering dimensions 
of the operator and its angular momentum. Indeed, all listed operators belong to 
irreducible representations of the angular momentum. 

On the lattice, the discretization of the covariant derivative can be done in a 
standard way: 


V/.'0(x) =-[Uf,{x)'il!{x + ajx) - 'iIj{x)] 

CL 

\7lip{x) =-[ijj{x) - Uf,{x - afi)~^'ip{x - afi)] 

CL 


( 17 ) 


but in general, owing to the lower (hypercubic) symmetry of the Euclidean lattice 
action with respect to that of the continuum ( all 4-d rotations), the identification 
of a given irreducible representation may require some particular combinations of 
operators. This classification has been discussed for example in refs. JOT and |TT 


With the Schrodinger functional the symmetry is in general further reduced, from 
hypercubic to cubic, because of the fixed boundary conditions in time that mark 
this direction with respect to the others. In this paper, we will concentrate on the 
calculation of the first two moments to which we associate the following irreducible 
operators: 


1 - ^ X-i 

0\2 T>2] y '^ 

1 - ^ ^ A-f 

Ol 23 =gV’7[i C 2 -D 3 ] yi/'- 


( 18 ) 


where is a flavour matrix. These are a subset of the basis described in ref. JIT 
involving only spatial indices and are multiplicatively renormalizable. We define the 
matrix elements of the first two moments by the observables: 
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(19) 


= -a® “>{i<^(x) 7 [i D21 jT®i/i(x)C(y)riT®C(z)) 

y,z 

/oi23(a^o) =f3{xo) = D 2 D 3 ] ^r^^(x)C(y)^rr3C(z)) 

y.z 

where the contraction of the classical helds is non-vanishing if the matrix F satishes: 

= -P+(-) and p is the momentum of the classical held sitting on the boundary. 
They can be seen as the operator matrix elements between the vacuum and “p”-like 
classical state sitting at the T = 0 boundary. 

We take the limit of massless quarks, which in the numerical simulations can 
be monitored via axial Ward identities. In the SF framework it is possible to work 
at zero physical quark mass because a natural infrared cutoff to the Dirac operator 
eigenmodes is provided by the time extent of the lattice. This choice simplihes 
the recursive procedure at hnite volume that leads to the reconstruction of the 
continuum non-perturbative renormalization constant of the operators. 

The matrix element of the operator for the hrst moment involves two directions 
and three for the second moment. These directions must be provided by external 
vectors: we have chosen to obtain one of them from the contraction matrix T, i.e. 
from the polarization of the vector classical state: 

r = 7 fc, A; = 1 , 2,3 (20) 

and the remaining ones from the momentum p of the classical Fermi held at the 
boundary. This choice gives for the tree level a non-vanishing matrix element in the 
massless quark limit, where we evaluate our correlations. The tree-level correlation 
can be easily calculated and reads, for the hrst moment: 

/f (a;o) = ^ [i-tpo) (M_(p)e-2-(p)"o - M+(p)e-2-(p)(2^-o))] (21) 

where 


Pf, = (2/a) sm{ap^/2), 

p^ = (1/a) sin(ap/,) 

(22) 

M{p) 

= m + 

^ -2 

M± = M{p) ± ipo 

(23) 

sinh 

[ 5 ^( 9 )] 

a ( 

+ {m+ ^aq^)^ | ^ 

(24) 

^ 2\ 

1 a{m + |aq2) / 
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R{p) = M{p) {1 - -ipo{l + . (25) 

We have chosen for convention F = 72 and p = (pi, 0,0). In the continuum 
and in the massless limit this expression reduces to: 

= (1 . (26) 

Notice that the symmetric choice F = 72^2 + 7 ipl and p = (pl,p2,0) would 
lead to: 


/IVo) (4e-"''-'7 

^ PiP 2 M{p) N /^^-2u>(p)xo _|_ 2 g- 2 <^(p)( 2 T-xo)'\ 

R{piY ^ 


(27) 


PiP2M{p){-ipo)N ^^^_ 2 a;(p)^„ 

r{pY ^ 


2e-2^(p)(2'r-a;o)^ 


which in the massless continuum limit reduces to 


f2°\xo) 


pIp2N /4g-2|p|T^ , PipIN 

(l + e-2|p|r)2 V ^ ^ (l + e-2|pl'r)2 


(4e-2|p|r) 


(28) 


In the massless continuum limit, only the contribution proportional to exp(—2cc;T) 
survives, while at hnite lattice spacing and for xq smaller than T the dominant term 
is exp(—2cc;a;o), i.e. in this case the signal is dominated by lattice artefacts. Another 
advantage of using the non symmetric contraction is the economy in the number of 
non-zero momentum components. 

The quantization of momenta on a hnite lattice is one of the major sources of 
lattice artefacts. Indeed, with ordinary periodic boundary conditions the momentum 
cannot be chosen, in each direction, smaller than a minimum value proportional to 
the inverse lattice size. The presence of more independent momentum components 
increases the absolute value of the total momentum and therefore of the lattice 
artefacts associated to its size. In this respect, we will see that a different choice 
of boundary conditions can provide a “hnite-size” momentum not submitted to 
quantization. 

For the second moment, we choose F = 73 and p = (pl,p2,0). The tree-level 
expression in this case is just proportional to the one for the hrst moment: 


ff\xo) = [(*Po) {M_{p)e - M+(p)e 2‘^(p)( 2T a;o)^j ^ ip 2 f^^\xo) 

(29) 
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The observable is taken as the ratio of the correlation at one loop normalized by 
its tree-level expression and mnst be expnrgated from the renormalization constant 
of the classical bonndary helds (. Following ref. []^, this is achieved through the 


normalization by the quantity called /i, also normalized by its tree-level expression. 

Finally, some care should be taken so as to ensure the massless quark limit at 
order The breaking of chiral simmetry of the Wilson action entails a non-zero 
shift of the quark mass from the naive value at order 


= — 


2 r d^q I 1 

3 J_„ 2tt^ \ 


3 + 


- W) 


e 




( 30 ) 


Including the mass shift simply amounts to a further contribution coming from 
the derivative of the tree-level expression calculated at the zero naive quark mass 
value times the coefficient of the order g"^ mass shift. 


4 The results of the calculation 

The expression for the observable at one loop in the continuum can be parametrized 
as: 


Z(pL, xo/L, a/L) = 1 + g'^Z^^\a/L) 

ryil)/ /r\ 7 I 1 / /r\ I CklTl(a/ L) ( 31 ) 

With \a/L) = &o + Co ln(a/L) -f 2_^ o- - jj. -, 

k=l 

where the scaling variables pL and Xq/L are kept constant when the number of lattice 
points is sent to inhnity. Indeed, beyond the free field case, this limit does not exist 
because of the ultraviolet divergence manifested through the logarithmic term that 
breaks the free-held theory scale invariance. In the following, the momentum will 
be set to its minimum value, i.e. pL = 2 ti. 

The operator needs renormalization to be hnite in the continuum limit: for 
example one can £x the operator matrix element equal to its tree-level value at 
p = 1/L. At one loop one gets: 


0 “^ =(1 - g^Z'')(a,,))0‘“”(a) 

=(1 - g"Z(‘>(a,7))(l + j"Z<‘>(o/L))0"'". 

The coefficients Cq are the interesting ones in the continuum limit, but at 
hnite lattice spacing, i.e. with a hnite number of lattice points, the correlation also 
contains corrections that decrease with the inverse powers of the number of points 
N. Our results are for the ordinary Wilson action without any improvement, and 
lattice artefacts start at order 1/N for both coefficients. 
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General formulae for calculating the perturbative expansion within the SF can 
be found in refs. |^, where the correlation in eqs. (^to H) is calculated as: 


[^(a^)C(y)]F = 
[^(a;)C'(y)]F = 


Sjjcljx) 

V(y) 

S'ii’cijx) 

<5p'(y) 

dp(x) 




5p'(y) 


(33) 


and in turn 


^ci(x) = 0 ^^^ {S(x,p)l/(p - ad,0) V(y)lyo=a 
y 

+ S(a:,!/)[/(!/, 0)p'(y)|„.T-a} 


( 34 ) 


We have done the calculation independently by using the notion of 't/^ci and the 
explicit expressions of ref. 0 or by using directly the correlations given by eqs. ^ 
to ^ and the perturbative expansion of the Dirac propagator. In the latter case we 
obtain also diagrams where the gluon field emerges from the time links connecting 
the Dirac propagator to the zero time slice. In both cases the calculation was done 
in a time-momentum representation. The algebraic manipulations leading to the 
finite sums performed numerically were done with the algebraic program “FORM”. 

The final Fortran result was obtained with double-precision accuracy. We have 
run all even lattice sizes ranging from 8 to 36 in the case of the first moment and to 
40 for the second moment. The latter choice was forced by the slower approach to 
the continuum for the second moment. The £t to the N dependence of our results 
has been made by the expression 


Z(^\N) = Bo + Coln(N) ^ )With N = L/a, (35) 

k=l,2 

excluding higher-order terms. The stability of the £t was checked against a change 
in the value of the smallest volume included. After checking the agreement to the 
per cent level with the expected value of the coefficient of the logarithm Cq, we fix 
it to its theoretical value and determine the remaining parameters. As an example, 
we report in Table 1 the results for the first moment and for the correlations at 
Xo = T/2 of the constant Bq. 









Range 

Parameters 

Constant 


24-34 

5 

0.0971 

10 -^ 

22-34 

5 

0.0968 

10 -^ 

20-34 

5 

0.0964 

10 -^ 


Table 1: Values of the constant of the correlation function of the first moment with xq = L/2. 


Range 

Eliminated parameters 

Constant 


24-30 

ln(Af)/Af, 1/Af, l/N'^ 

0.0967 

1.5 

24-30 

\n{N)/N, l/N 

0.0967 

o 

1 

to 


Table 2: Values of the constant of the correlation function of the first moment with xq = 
after elimination of some subleading contributions. 


The values are shown for different fitting intervals and, in the last column, we 
report the values with a reference error of 10“® on the theoretical points . An 
alternative determination can be achieved by making suitable combinations of the 
data where some of the dominant 1/L and 1/L^ corrections are eliminated. This 
leads to the results of Table 2, where we indicate the type of fitting coefficients 
eliminated from the general eq. and the fitting interval used. 

Our final results for the constants, for the first and second moment, are in 

Table 

The approach to the continuum values is very slow, as can be seen from figures 
|l| and H where we compare the ratio of perturbative finite lattice coefficients with 
their continuum values for the two moments and the two xq cases. 

In the case of the second moment the slope of the approach is lower, indicat¬ 
ing that momentum quantization is an important source of lattice artefacts. The 
continuum is approached from above in all cases, similarly to what happens for the 
proton matrix elements, by comparing the standard Wilson action with a partially 
non-perturbative 0{a) improvement |^. 

The use of a physical momentum is not the only way to obtain non-vanishing 
matrix elements for the operators that we have considered. Indeed, the choice of 


Moment 

Definition 

Constant 

First 

Xo = T/4 

Bq = 0.084(1) 

First 

Xo = T/2 

Bo = 0.0967(5) 

Second 

Xo = T/4 

Bo = 0.010(5) 

Second 

Xo = T/2 

Bo = 0.008(2) 


Table 3: Values of the constants of the correlation functions A and /a with the operator insertion 
at xo = L12 and xg = LjA. 
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Figure 1: The ratio of perturbative finite lattice renormalization constant with their continuum 
value for the first moment in the two cases of Xo = L/4(o) and xo = L/2{»). 


periodic boundary conditions for fermions different from the usual ones 


P 


to: 


'ip{x) = 'ip{x + Lk) 


(36) 


il){x) = + Lk) (37) 

where 9 is an arbitrary constant phase, introduces a “hnite-size momentum 9" 
probed by the operator only when it involves helds on the two sides of the boundary 
that can play the same role as the ordinary momentum. One can also distribute the 
phase to all lattice points by an Abelian transformation on the Fermi helds: 

'ip{x) = e^i^ip{x), (38) 

which changes the form of the lattice derivative from 
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2 



Figure 2: The ratio of perturbative finite lattice renormalization constants with their continuum 
value for the second moment in the two cases of xq = L/4(o) and xq = Lj^i*). 


=-[U^{x)'^{x + afi) - 'i/j{x)] 

1 _ . 

X/l'ipix) =-['0(a;) - Uf,{x - afi)~^'ijj{x - afi)] 

CL 

to: 

\7^,'ip{x) =-[\^U^{x)'ip{x + ap,) - p{x)] 

“ (40) 

=-[ij{x) - Xp^U^ix - ap)~^ij{x - ap)], 

CL 

where 

6o = 0, -TT < Ok < TT. (41) 
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In this case the operator feels a “hnite-size momentum 9/L" aX all lattice points. 
By taking the average of the operator insertion over all points, the same result is 
recovered. 

We have repeated the perturbative calculation at zero momentum and non¬ 
zero “hnite size momentum” 0 = 0.1/L, to be compared with the minimal lattice 
momentum = Stt/L. 

The lattice results converge much faster to their continuum limit, confirming 
the momentum quantization as a major source of lattice artefacts. Figures ^ and ^ 
show the ratio of lattice over continuum results also for this case. 



Figure 3: The ratio of perturbative finite lattice renormalization constants with their continuum 
value for the first moment in the two cases of xo = L/4(o) and xq = Lj^i*) using “finite size 
momentum”. 

In this case also a global fit without elimination of subleading terms provides 
a very good estimate of the value of the constant. Our final best estimates are given 
in Table 

The result for the constant of the continuum limit of the quantity /i that 
subtracts the renormalization constants of the C fields from that of the operator is 
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o x_0 = L/4 
• x_0 = L/2 




1 







0.5 


Q I _^^^^^^_ I_I _^^^^^_ I 

12 3 4 

Log(L) 


Figure 4: The ratio of perturbative finite lattice renormalization constants with their continuum 
value for the second moment in the two cases of xq = L/4(o) and xo = L/2[*) using “finite size 
momentum”. 


Bq = -0.3590(1). (42) 

It has been obtained from fits to the data made available to us by P.Weisz. 

Table ^ contains a summary of the constants Bq for the operators, after remov¬ 
ing the external legs renormalization, in the various cases that we have discussed. 

We have shown by an explicit one-loop calculation the feasibility of a lat¬ 
tice evaluation of the renormalization constants of the first two non-singlet twist-2 
operators within the Schrodinger functional scheme. The definition introduced in 
this paper is suitable for a finite-volume recursive scheme that can provide a full 
numerical reconstruction of the scale dependence of the renormalization constants, 
even in the non-perturbative regime. By comparing lattice and continuum pertur¬ 
bative estimates, we can conclude that lattice momentum quantization in a hnite 
volume is one of the most important sources of lattice artefacts. The modihcation of 
fermion boundary conditions already introduced in the framework of the Symanzik 
improvement programme carried out by the ALPHA collaboration, introduces of a 
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Moment 

Definition 

Constant 

First 

Xq = LjA 

Bq = -0.0577(1) 

First 

O 

to 

Bq = -0.4129(1) 

Second 

Xq = LjA 

Bq = -0.2403(5) 

Second 

Xq = L/2 

So = -0.747(2) 


Table 4: Values of the constants of the correlation functions /2 and /a in the two definitions 
with the “finite size momentum”. 


Moment 

Definition 

Constant 

First 

Xq = L/A{p ^ 0) 

Bq = 0.2635(10) 

First 

Xq = L/2{p ^ 0) 

Bq = 0.2762(5) 

Second 

Xq = L/A{p ^ 0) 

Bq = 0.1875(20) 

Second 

Xq = L/2{p ^ 0) 

Bq = 0.1895(50) 

First 

Xo = L/A{e ^ 0) 

Bq = 0.12180(15) 

First 

Xq = L/2{e ^ 0) 

Bq = -0.23340(15) 

Second 

Xo = L/A{e ^ 0) 

Bq = -0.06080(55) 

Second 

Xo = L/2{9 ^ 0) 

So = -0.5675(20) 


Table 5: Values of the constants of the operators renormalized in the two definitions with the 
real momentum p or the “finite size momentum” 0 different from zero. 


“finite-size momentum”, which escapes the quantization rule and turns out to sub¬ 
stantially improve the approach to the continuum of our lattice results. Numerical 
studies with standard lattice momenta and with finite-size momenta are under way 

fPH. 
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